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Abstract 


This work addresses the problem of optimal positioning of a set of measured points 
with respect to an ideal design surface. Localization refers to the process of determining the 
rigid body translations and rotations which must be performed on the set of points to move 
those points into closest correspondence with the design surface. In unconstrained 
localization all points have equal effect on the determination of the rigid body 
transformation, while constrained localization allows a subset of the points to have stronger 
influence on the transformation. 


The measured points 1n the context of this work refer to physical points in space that 
are obtained by direct measurement of a manufactured marine propeller blade. The ideal 
design surface refers to a surface descnption of the propeller blade provided by the blade 
designer. Given that the measured blade is manufactured from the design surface 
description, it is the task of localization to determine an optimal positioning that will bring 
the measured points of the manufactured surface as close as possible to the design surface. If 
the manufactured blade is repositioned in space according to the prescription of the 
localization transformation, it will have the closest possible correspondence to the original 
design. Direct benefits to the manufacturer may result from less wasted material in initial 
castings and better ability to program postcasting work through optimal positioning of the 
workpiece. 


The constrained and unconstrained localization method is developed from a theoretical 
basis. Applications of the localization method are investigated with examples of propeller 
designs and inspection data obtained from blades that were manufactured from those designs. 
Experimental results demonstrate the capabilities of the method and its applicability to 
automated inspection. 
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Chapter 1 


INTRODUCTION 


A fundamental problem in manufacturing is the need to determine if a manufactured 
piece meets the requirements of the original design description from which it was made. The 
evaluation of positional tolerances to ensure that a manufactured item is an acceptable 
rendering of the original design is a basic element of manufacturing inspection. 

In few areas of manufacturing inspection is the need for precise inspection and 
evaluation of tolerances more clearly demonstrated than in the area of marine propellers. 
The manufactured item is an exceedingly complex sculptured surface which must be 
produced with extremely high fidelity to the original design. Very strict positional tolerances 
must be achieved in a difficult manufacturing process to prevent severe compromise of the 
performance of the propeller. 

The inspection of marine propellers has traditionally involved highly skilled 
technicians checking the surface of a manufactured propeller with numerous mechanical 
gages. The gages are cut to specified dimensions by the direction of a manufacturing 
engineer who interprets the specifications of the propeller designer. Although rigid 
guidelines are provided for placement of the gages on the blade, errors in measurement can 
result from decisions by the technician regarding "fit" or alignment of the gage on the 
manufactured blade. Moreover, the direct gage measurements only evaluate the blade at the 
local site where the measurement is made. A completely satisfactory method for evaluating 


global compliance with specified tolerances has not been available. Often expensive and 
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inefficient rework of propellers has been necessary because it was not possible to quickly and 
confidently ascertain whether the manufactured product would satisfy the requirements of the 
designer. 

The recent development of automated methods of inspection using coordinate 
measuring machines (CMM) and laser-based measuring devices has made it possible to 
obtain voluminous quantities of highly accurate spatial measurements of manufactured 
propellers. These robotic devices have provided a reliable source of measurement data, but 
methods for best using that data are still being developed. 

This thesis will address an aspect of the question of how to best utilize measured data 
from manufactured propellers for the automated inspection of those propellers. It will deal 
with the problem of optimal positioning of a set of measured points from a manufactured 
propeller blade relative to the design surface from which the blade was manufactured. The 
problem investigated in this thesis may be simply stated as follows: 


Given a set of measured data points from a manufactured surface, 
determine the rigid body translations and rotations which must be applied to the 
set of measured data points to bring those points into closest correspondence 


with the design surface from which the measured surface was manufactured. 


If all measured points contribute equally to the determination of the set of 
rigid body transformations, then the problem is one of unconstrained 
If some measured points have greater effect on the determination of the set 


of rigid body transformations than other points, then the problem is one of 
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The thesis will develop the investigation of this problem using the following structure. 

Chapter 2 wil present a review of current literature relevant to the problem of 
localization. Particular emphasis will be given to literature which is directly related to the 
development of the localization algorithms presented in this thesis. 

Chapter 3 will discuss the theoretical basis for the unconstrained localization algorithm. 
It will describe the optimization problem and the procedures that are followed to obtain a 
rigid body transformation matrix which is a solution to the unconstrained localization 
problem. 

Chapter 4 will develop the theoretical basis for the constrained localization algorithm. 
It will contrast the differences between constrained and unconstrained localization and will 
describe the procedures involved in the solution of the constrained problem. 

Chapter 5 will present experimental results which demonstrate the applicability to the 
problem of localizing manufactured marine propellers. One application involves the 
localization of a simple fan blade, while another application involves the localization of a 
complex marine propeller blade. In each case results are presented to demonstrate the 
usefulness of the constrained and unconstrained localization methods. 

Chapter 6 summarizes the results of the investigation and presents ideas for additional 
work. 

Appendix A will demonstrate the orthogonality of the rotational transformation matrix 
that is used in this investigation. 

Appendix B will derive the Jacobian matrix for the objective function used in the 


unconstrained localization algorithm. 





Appendix C will use orthogonality to develop a method for determining the parameters 
in a design surface of the projection of a point onto that surface to a very high degree of 
accuracy. The method finds direct application in improving the accuracy of the calculation 


of the minimum distance from a point to a parametric surface. 
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Chapter 2 


LITERATURE REVIEW 


There has been much interest in the problem of localization in recent years and the 
current literature reveals some aspects of methods that have been used in developing 
solutions to the problem. This chapter will review some of the literature that is relevant to 
the localization problem and the work of this investigation. The intent is to provide some 
pertinent background information which will give the reader a broader perspective 
concerning this particular work. 

Localization of surfaces was accomplished by Thome using a least squares matching of 
associated boundary edges [Thome 85]. Gunnarsson developed a localization method 
between a set of points and parametric surfaces by dynamically faceting the design surface 
and finding a rigid body transformation matrix which minimizes the sum of squared 
distances from the data points to associated planar faces [Gunnarsson 87a], [Gunnarsson 
87b]. This formulation required the solution of a constrained minimization problem with 12 
unknowns, which are the elements of the rotation matrix with three translations, and 6 
constraints which are the necessary relations between elements of the Euclidean rotation 
matrix. 

One question that is intrinsic to the formulation of the localization problem is that of 
selection of an appropriate norm to use in distance minimization. Bourdet and Clement 
present an analysis of this problem in [Bourdet 88]. They show, through the use of a 
small-displacement screw linearization model, that the L,, (or minimax) norm reduces form 


error when compared to the L, (or least-squares) norm by approximately 15 per cent for a 
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small number of points. The improvement however disappears when the number of points is 
greater than twenty. In contrast, the L, norm is superior to the L,, norm in its usefulness for 
statistically detecting aberrant points in a particular data set. Moreover, the L., norm was 
shown to require much greater computational tume than the L, norm. The general guidelines 
of the paper can be summarized by saying that for point sets of fewer than 12 points the L. 
norm should be used, while the L, norm should be used for point sets of greater than 12 
points. An interesting algorithm which allows for computation of both the L, and the L, 
norms at the same time is presented in [Goch 90]. 

Some of the theoretical framework upon which this work is based is presented by Alt, 
Mehlhom, Wagener and Welz] in [Alt 88]. The writers of the paper demonstrate some 
algorithms which pertain to the problem of mapping congruent objects A, B in KR” from one 
to the other using geometric transformations involving rotations, translation, reflections, and 
stretching. Although the algorithms developed in this thesis do not directly address reflection 
and stretching, the mapping problem 1s essentially the same one. A significant finding in 
[Alt 88] is that exact congruencies are not possible using measured data. Even small 
perturbations in the measurements will destroy conguencies between the two geometries. 
They therefore establish the approximate congruence problem with tolerance €, where the 
maximum distance between two corresponding sets of points is less than or equal to €. This 
problem is shown to have a computational upper bound of O(n’ log n) for either L, or L.. 
norms, where 7 is the number of points, when the sets of points are assumed to be known and 
a congruency relationship exists between the two point sets. The writers of the paper deal 
only with point sets; they do not discuss the particular problem of mapping approximately 


congruent point sets to a surface. 
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Imai, Sumino and Imai continue work on the problem of mapping point sets in [Imai 
88]. They develop a minimization algorithm which greatly improves the theoretical 
efficiency of the methods of Alt, et al. They show that, by using their algorithm, the L, norm 
problem has a computational upper bound of O(n’) and the L. norm problem has a 
computational upper bound of O( n+n log n) in the worst case, where 7 is again the number 
of points. 

An essential element of the localization algorithms presented in this thesis is the use of 
orthogonal projection methods of curves onto surfaces to locate nearest points in the 
parametric space of a design surface. These methods were developed originally by Pegna 
and Wolter in [Pegna 90]. In this paper a space curve is mapped onto a surface by tracing a 
surface curve whose points are connected to the space curve by surface normals. The tracing 
is achieved by solving a tensonal differential equation in the parametric space of the surface. 

Another concept that is important in this work was developed by Kriezis, Patrikalakis 
and Wolter in [Kniezis 90], [Knezis 91] for use in the solution of surface intersection 
problems. It is the oriented distance function which is defined as the inner product between a 
vector from a given point to its nearest point on a surface, and the unit normal vector of the 
surface at that point. The method is used in formulating the constrained localization problem 
of this investigation under the assumption that the two vectors are collinear. 

The fundamental basis for the work of this thesis applied to the unconstrained 
localization problem derives from earlier work presented by Bardis and Patrikalakis in 
[Patrikalakis 90]. In this work positional tolerances were represented in terms of a ball offset 
tolerance region around an ideal rational spline surface (design surface). The tolerance 
region bounding surfaces were approximated by rational B-splines. The manufactured 


surface (target surface), known either in terms of a lattice of measured points or as a 
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simulated surface from numerical or analytical predictions, was then optimally positioned 
with respect to the ideal design surface by minimizing the L, distance norm. After 
localization, the target surface could be tested for intersection with the boundaries of the 
tolerance region. If the target surface is found to be entirely within this tolerance region, then 
the agreement would be considered satisfactory. These unconstrained localization techniques 
were further developed to improve their computational efficiency, and they are presented in a 


paper coauthored by this writer [Bardis 91]. 
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Chapter 3 


UNCONSTRAINED LOCALIZATION 


3.1 Introduction 


Localization can be defined as the problem of determining the optimal positioning 
of a set of measured points relative to a design surface. If all measured data points have 
equal effect on the determination of this optimal positioning, then the localization 
problem can be defined as a problem of unconstrained localization. 

The localization problem can be formulated as an optimal parameter estimation 
problem involving six parameters. Those six parameters are the three translations and 
three Euler angles which correspond to a general three-dimensional translation and 
rotation of a rigid body in space. The problem can be formulated as an unconstrained 
minimization, where the objective function of the minimization is the sum of squared 
minimum distances of a set of measured points from a design surface. In this context the 
measured points represent physical points in space that are determined by direct 
measurement of a manufactured surface. The design surface refers to the underlying 
design description that is used to produce the manufactured surface. The unconstrained 
localization problem is then the problem of estimating the six parameters of a rigid body 
transformation which will bring the set of equally weighted measured points into closest 


correspondence with the design surface. 
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3.2 Problem Formulation 


Consider a parametric surface P(u,v), which will be called the design surface, 


representing the desired design geometry. Consider also a set of m points R,, 
1<ism, R, € &’, which will be called measured points. Finally, consider another set 
of m points Q;, 1 Si <™m, in the design surface P(u,v) which are the nearest points to 
each measured point R;. It is assumed that the nearest points Q; are unique and are the 
grthogonal projections of the points R; onto P(u,v). The points Q; will be subsequently 
called projections. 

The assumption that the nearest points Q; are unique and also orthogonal 
projections of the points R; onto P(u,v) is reasonable for the inspection problem that is 
being investigated. The measured points of a manufactured propeller blade can be 
expected to be very close to the surface of the propeller and relatively far from the medial 
axis of the blade. It can be shown that under some reasonable conditions, the nearest 
point will be unique [Kriezis 91]. Furthermore, the only places on the blade where the 
orthogonal projection condition could be expected to not be valid would be at the trailing 
edge, root edge or tip of the blade, and these locations are ones where measurements 
cannot be readily made anyway because of the edge discontinuity. ' 

The objective function for minimization is the sum of squared minimum distances 
of each point from the design surface. Each measured point R; has a nearest point Q; on 


the design surface. The minimum distance from each measured point to the design 


'Even though the assumption of the existence of a unique nearest point which is also an 
orthogonal projection of a measured point to the design surface is not unrealistic for this 
problem, its importance cannot be overemphasized. It forms a foundation for much of the 
theoretical development of this thesis. 
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surface can be simply defined as the Euclidean distance between the points R; and Q,. If 
the distance between two points is denoted by d(P,, P,) then the minimum distance from a 


measured point R; to the design surface P(u,v) is defined as 


d(R;,Q;) = | R;-Q;| = min, ,d [R;, P(u, v)] Cy) 


The squared minimum distance is 


(R,; - Q;) eA) SS Q,) = [d(R,, Q,)} (3-2) 
and the objective function becomes 


OF = ¥[d(R,,Q)F (3-3) 
The parametric surface is assumed to be an untrimmed rational B-spline (NURBS) 
patch of orders M in u, and N in v, (OS u,v <1). The surface is further assumed to have 
continuous first order partial derivatives (C' continuity). The design surface patch is then 


of the form 
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3 a h;B; m(u)B;.w(v) 


imO j 
PO sla 10 =) S721) ate the vertices of the associated control polyhedron and 
h;,(0Sisgm-—1, O0OSjSn-1) are positive weights. B;,(u) and B;v) are the B-spline 
basis functions over open knot vectors [Gordon 74] with variable knot spacing of 
Oho ER een Ser me (3-5) 
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The goal of the localization problem is to obtain an optimal rigid body 
transformation for operation on the set of measured points, R,, so that those points will 
correspond as closely as possible to the design surface, P(u,v). 

It is possible to generalize the objective function to accomodate multiple untrimmed 
patches using analogs of the equations for a single patch. Consider n patches such that 


P(u,v), 1SjSn, represents the jth patch, and R,;, 1Sism, lSjsn, Rie ®’, 


uy? 
represents the ith measured point associated with the jth patch. If Q;; now represents the 
minimum distance orthogonal projection of R; onto the design surface P;(u,v), then the 


squared minimum distance for multiple patches becomes 


(R; ~ Q;) ; (Ri, = Q;)) = [a(R;,Q,)} (3-/) 


and the objective function for multiple patches becomes 


OF = Y(dR;,Q)F (3-8) 


3.3 Localization Transformation 


The localization procedure that is used to determine the optimum rigid body 
transformation for the measured points consists of minimizing the objective function of 
(3-3) by calculating values for the three rotations and three translations which are the six 
parameters of a rigid body transformation. 

The three rotational parameters of the mgid body transformation are given by the 


angles yw, 8 and @ (Euler angles) which represent rotations about the x, y and z axes 
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respectively. The three translational parameters are represented as elements of a 
translation vector t, having components ¢,, ¢, and ¢, for translations along each respective 
axis. 

If the set of measured points of the localization problem is operated upon by 
successive rotations followed by successive translations, then a new set of points, r,, can 
be obtained from the original measured points by the operation 

r= [C]R,+t (3-9) 

In equation (3-6) the matrix [C] is defined as the rotational transformation 
matrix, which is obtained by multiplying the matrices associated with rotations about 
each coordinate axis. If rotations are ordered in the sequence rotation by angle @ about 
the z-axis, followed by rotation by angle 8 about the y-axis, followed finally by rotation 


by angle y about the x-axis, then the matrix [C] 1s given by 


l 0 0 cos8 QO sin@\{cosd —sind 
[C] = |0 cosy -siny 0 ] 0 sind cosd QO (3-10) 
O siny cosy /\-—sin8 O cos®@ 0 0 ] 
or 
cos 8cos —cos 8sind sin 8 
[(C] = | cosysind+sinO9sinycosd cosycosd—sinOsindsiny  —sinycos® (11) 


sinysind—sinO@cosdcosy sinycosd+sinOsindcosy  cosOcosy 
The matrix [C] is an orthogonal matrix, having the property [C][C]’ =[C]' [C]=4,, 


where /, is the identity matrix of dimension 3.’ 


* A demonstration showing that [C] is orthogonal can be found in Appendix A. 
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3.4 Localization Algorithm 


The localization algorithm is the unconstrained minimization of the objective 
function given as equation (3-3). The process involves iterative operations on the set of 
measured points R,, producing new sets of transformed points according to equation 
(3-9). If r; is a transformed point given by equation (3-9), and q,; is the minimum distance 
orthogonal projection of the transformed point onto the design surface, then the minimum 
distance squared from r; to P(u,v), [d(r;, q)y, is defined by the following equation which 
follows directly from equation (3-2): 

[d(r,,q,)P = (r,-4,) * (t;- 4) a) 
Using the transformed points r; and the corresponding projections q; to calculate the 
squared minimum distance at each iteration step, the iterative process continues until a 
minimum is reached. This minimization can be simply stated as the following: 


Determine 0, 9, y,¢,,¢,,¢,, such that 


xy? 


m (3-13) 
OF (, 9, W,t,5t,56,) = 2{4 (r;,q;)) 
iS minimized. 
An alternative objective function would be the minimax (or Tschebyscheff norm), 


L.. , where 


OF (4, 3) Leet) = max; | Sup [d(r,;, q;)] — Inf [d(r;, q;) Ill (3-14) 
is minimized instead of equation (3-13). Such a norm allows calculation of the 


parameters minimizing the maximum of the minimum distance of all measured points 


from the design surface. However an objective function of this type is much harder to 
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implement for a large number of points, and [Bourdet 88] has shown that no significant 
improvement over the L, norm results if the number of measured points is greater than 
about 20. 

Determination of the minimum distance from a point to a parametric surface using 
(3-1) involves the calculation of a minimum with respect to the design surface parameters 
u and v. A modified Newton algorithm implemented in the Numerical Algorithm Group 
routine EO4KCF, [Gul 76], [NAG 89], is used with an initial guess of the minimum 
distance position provided by using the u, v parameters of an orthogonal projection onto 
Py (u,v) of B-spline curves that are fitted through selected sets of data points. This 
method uses the orthogonal projection techniques of curves on surfaces developed in 
[Pegna 90]. It has resulted in very good computational efficiency for two reasons: 


1. An exhaustive search of Py y(u,v) to compute minimum distance needs to be 


performed only once for the end point of each B-spline curve to find its 
minimum distance orthogonal projection. 

2. The B-spline approximation for the measured data can be very rough. 

The modified Newton algorithm used for calculation of the distance of each point to 
the design surface behaves well if r; is close to Py (u,v). For the small rotations and 
translations that are performed on each point using (3-9) it is unlikely that several local 
minima of d(r,;, Psy y(u,v)) will interfere with the process. 


The unknown quantities 6,9, y,7,,f,,f, which render the objective function (3-2) 


minimum are computed by applying a quasi-Newton algorithm implemented in 


Numerical Algorithms Group routine EO4JAF. Estimates of the Jacobian and Hessian of 


ie 





OF ($, 9, Y,t,,t,,t,) are used to generate a series of feasible sixtuples converging to a 
minimum, [Gill 74], [NAG 89].° The initial guess of the unknown quantities is chosen as 
the zero vector. 


After one calculation of rotations and translations, new points r; are created using 
equation (3-9). The process continues until a minimum for the objective function is 
achieved within a tolerance level specified by the user. At the conclusion of the process 
the user will be provided with an optimal prescription for the three rotations and three 
translations which should be applied to the measured points to bring them into closest 


correspondence with the design surface. 


* The Jacobian of the objective function is estimated in the NAG routine EO4JAF by 


difference quotients. Direct calculation of the Jacobian might however improve the 
performance and accuracy of the quasi-Newton algorithm if implemented in another routine. 
The derivation of the Jacobian for the objective function is therefore provided in Appendix 
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Chapter 4 


CONSTRAINED LOCALIZATION 


4.1 Introduction 


Unconstrained localization addresses the problem of determining the optimal 
positioning of a set of measured points relative to a design surface when each measured 
point has an equal effect on the determination of positioning. Constrained localization, 
on the other hand, involves the problem of determining a feasible, but not necessarily 
optimal, positioning of a set of measured points relative to a design surface when subsets 
of the measured points can have unequal effects on the determination of positioning. 

The unconstrained localization problem seeks to minimize one global objective 
function so that measured points are all collectively brought as close to the design surface 
as possible. In this problem each point contributes with the same weight to the 
minimization of the objective function. In contrast, the constrained localization problem 
Starts with the rotation and translation produced by the minimized objective function of 
the unconstrained localization problem and determines a rigid body transformation which 
will allow the measured points to satisfy a set of nonlinear constraints. The constrained 
localization problem does not minimize an objective function, but rather uses 
minimization techniques to find a feasible transformation that will satisfy the constraints 
imposed by a set of constraint functions. Satisfying the constraint functions has the effect 


of changing the importance of each measured point. 
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The need for constrained localization of marine propellers is evident because of the 
tighter positional tolerances that are required near the leading edge of a blade. A method 
is required which will provide for greater influence on the localization by measured 
points that are close to the leading edge. 

The typical design of a propeller blade provides a convenient way to construct the 
constrained localization problem. Normally the blade is described by the designer in 
terms of a NURBS surface with two sets of isoparameter lines running approximately 
parallel and perpendicular to the leading edge of the blade. In this manner the 
isoparameter lines form a grid in the spanwise and chordwise directions of the blade. 
Two spanwise isoparameter lines can be assigned to bound the leading edge region of the 
blade. Constraints on the localization can then be imposed based upon whether or not 
parametric values of the minimum distance orthogonal projections of measured points on 
the design surface lie inside or outside of the boundaries defined by the two isoparameter 
lines. If the points within the two isoparameter lines in the parametric space of the 
surface are nearest points to the measured points, then measured points near the leading 
edge of the blade can have greater effect upon the localized positioning as required. This 
is accomplished by associating each measured point with its respective projection, and 
using the position of the projection in the parametric space of the design surface to 


determine the constrained function and consequent effect of the measured point upon the 


localization. 
The constrained localization method assumes that the region of the design surface 
that is associated with each measured point by its projection on the design surface doe 


not change during the localization, This assumption may not always be valid if the 


projections of measured points are very close to the isoparameter lines which bound the 
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leading edge region of the blade. For the small rotations and translations that are 
anticipated for this method, it is not expected that the projections of many measured 
points will change regions during the localization. Furthermore, even if such a 
phenomenon occurs the overall effects on the inspection problem are expected to be 


relatively small. 


4.2 Problem Formulation 


The constrained localization problem can be formulated as a nonlinear constrained 
optimization problem. The constraints are limits on distances of measured points from 
the design surface. They are determined by positions in the parametric space of the 
design surface of points which are nearest to measured points. In the formulation of this 
problem the constraints are inequality bounds on the constraint function which are 
allowed to have only two possible values.* If the points in the design surface lie within 
the isoparametric boundaries which define the leading edge of the blade, then the 
associated constraint on the optimization will have one value; if the points lie outside of 
this boundary, then the associated constraint will have another (larger) value. The 
problem is inhererently nonlinear because the rotation transformations which operate on 
the measured points to change their distances from the design surface are formed from 


combinations of transcendental functions in three independent variables. 


“Permitting only two values for the constraints reduces the complexity of the problem, but 
may make it more difficult to solve because of the forced discontinuity at the boundary. An 
alternative formulation would provide a transitional region for the constraint function to 
mitigate the transition boundary problem. 
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4.2.1 Global Objective Function 


It is assumed in the formulation of the constrained localization problem that the 
unconstrained localization problem has been previously solved using the techniques of 


Chapter 3. In this case the objective function given by 


OF = = [a,Ri} (4-1) 
has already been minimized. This means that at the beginning of the constrained 
localization an optimal global positioning to place the set of measured points at the 
minimum distance from the design surface has already been determined. Given this 
initial condition there is no need to further minimize a global objective function. On 
the other hand, constraint functions for each measured point must be determined, and 
a solution procedure must be performed to find an appropriate positioning which will 
satisfy the constraints. Since it 1s assumed that the constrained localization starts from 
a position of global optimization, it is expected that the constrained localization 
procedure will not produce an improvement in the global result obtained from the 
unconstrained localization. A satisfactory result will be to find a positioning of the 
measured points which satisfies the constraints. This result will not be necessarily 
unique. 

The assumptions associated with the global objective function for the 
constrained localization problem can be summarized in the following statements: 

(1) The constrained localization problem starts from a position of global 

optimization which is the solution to the unconstrained localization. 


(2) The constrained localization cannot improve the global result of the 
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unconstrained localization and therefore does not attempt a global optimization. 
(3) A satisfactory result is a solution, not necessarily unique, which satisfies the 


constraints. 


4.2.2 Constraint Function Definition 


The selection of an appropriate constraint function is fundamental to the 
formulation of the constrained localization problem. The constraint function must 
certainly be a distance measure, but careful definition of this measure may facilitate 


the solution of the problem. 


4.2.2.1 Squared Distance Function 


Consider again the parametric design surface, P(u,v), the set of m measured 


points, R;, and the set of m unique nearest points Q; on P(u,v), which are also 
assumed to be orthogonal projections. The minimum distance from a measured 


point R; to the design surface P(u,v) was defined as 

d(R;,Q;) = | R;-Q;| = min, _,d [R,, Plu, v)] (4-2) 
and the minimum distance squared was given as 

(R,-Q,)° (R,- Q,) = [4(R;, QP Gee) 
This function provides a measure of the proximity of a point to a surface and 
was chosen as the objective function for the unconstrained scree However 
the squared distance function is not the best one to use for every type of problem. 
For minimization problems involving very small changes in magnitude, the 


squared function introduces inaccuracy because all changes are squared. This 
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effect is manifested by reducing the number of decimal places that can be 
calculated with confidence by half. (From another point of view, double precision 
calculations effectively become single precision calculations.) Moreover the 
squared distance function as an unsigned function causes all sense of position of a 
measured point relative to the design surface normal vector orientation to be lost. 
The loss of positional sense of a measured point relative to the design surface is 
particularly undesirable in the context of using the process as an inspection tool. A 
manufacturing engineer evaluating the results of a localization operation would 
want to know whether he should cut or weld a manufactured blade. Therefore 
positional sense is essential. 

Good results were obtained for the unconstrained problem using the squared 
distance objective function, notwithstanding its previously mentioned 
shortcomings.” However it was judged at the beginning of the investigation of the 
constrained problem that the squared distance function might not provide 
sufficient numerical accuracy for solution of this more difficult nonlinear problem. 
An alternative distance function was desired specifically to enhance the numerical 


accuracy that could be expected from calculations. 


4.2.2.2 Oriented Distance Function 


A different method for determining the distance from a point to a surface 


using an oriented distance function was introduced in [Kriezis 90], [Kriezis 91]. 


> An expanded discussion of the results of the unconstrained localization is provided in 


Chapter 5. 
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This function does not use a squared quantity, and it retains the sense of relative 
position between a given point and a surface by using the normal of the nearest 
point in the surface to the given point. 

If n; is the unit normal vector of the design surface at the nearest point Q, on 
P(u,v), which 1s also an orthogonal projection of the measured point R; on P(u,v), 
then the oriented distance, d(R, Q), from R; to the design surface P(u,v) can be 
defined as 

d(R,,Q,) = n; + (R;-Q)) (4-4) 

This function will produce a very accurate measure of the distance from a 
point to a surface if the difference vector R;-—Q; can be calculated with high 
accuracy. A method for improving the accuracy of this calculation by exploiting 
the orthogonality of the difference vector to the design surface tangent plane is 
developed in Appendix C. 

The oriented distance function has a form similar to (4-2), but retains the 
positional sense of the unit normal vector n; at each projection on the design 


surface. 
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4.2.2.3 Constraint Assignment 


Using the oriented distance as the constraint function, the localization will 

have m constraints c;, one for each measured point R,, defined as 
c, = d(R,,Q,) (4-5) 
The constraints are assigned based upon the position of projections relative to the 
isoparameter lines which define the leading edge region of the blade using the 
following scheme. Assume the chordwise parametrization of the design surface in 
the parameter u with u;, u, the isoparameter lines at the leading edge boundaries. 
Let € be the value of the constraint for projections in the leading edge region and 6, 
the value of the constraint for other points. Typically, €«96. If u; is the u 
parameter in the design surface of the projection Q,, then constraints c;, can be 

assigned to each measured point according to 

—ESc,S+€ (4-6) 


fu, Su; Su, 


or —§ Sc, 548 (4-7) 


ifu; <u, OF u;>u, 


4.3 Constrained Localization Algorithm 


The constrained localization is the problem of determing the rotations and 
translations which must be performed on the set of measured points so that they will 
satisfy the required localization constraints. The problem can be summarized in the 


following problem statement. 
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For a set of m measured points R,, having nearest point orthogonal projections Q, on 
the parametric design surface P(u,v), where 1 Si Sm, determine the set of rigid body 
rotations and translations @, 9, Y, f,, 4, f, such that the oriented distances from R; to Q, 
satisfy the constraints c, as they are defined in (4-6) and (4-7) 

As in the case of the unconstrained localization problem, the measured points are 
operated upon by successive rotations followed by successive translations to produce a 
new set of points r;, defined by 

r = (CJR, +t (4-8) 
where t is again the translation vector with components ¢,, ¢,, t, and [C] is the rotational 


transformation matrix given by 


cos 8cos > —cos§@ sino sin 9 
[(C] = | cosysind+sinOsinycosd cosycosd—sinOsindsiny —sinycos@ (2) 
sinysind—sin@cosdcosy sinycosd+sinOsindcosy  cosOcosy 


The problem of determining the set of six parameters which will allow the set of 
measured points to satisfy the localization constraints is solved using the routine EO4UCF 
for nonlinear constrained optimization problems provided by the Numerical Algorithms 
Group (NAG). The routine uses an iterative sequential quadratic programming (SQP) 
algorithm in which the search direction is the solution of a quadratic programming (QP) 
problem, [Gill 86], [NAG 89]. 

The nonlinear constrained optimization routine estimates gradients of user-supplied 
functions with difference quotients unless the user can also supply those gradients. The 


latter situation produces a great improvement in computational accuracy and efficiency. 
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For this reason, part of the implementation of this algorithm involves providing symbolic 
gradients (or equivalently Jacobians), for each of the functions that are supplied to the 


NAG routine. 


4.3.1 Global Objective Function 


Although an objective function for the constrained localization is not minimized 
as previously explained, the structure of the nonlinear constrained optimization 
routine requires an objective function to be supplied. It is sufficient in this case to let 
the objective function be defined simply as a constant.° This is the simplest possible 
definition for an objective function and it allows the optimization routine to work 
accurately and rapidly. The objective function may be defined then as the following: 


Determine 6, 8, W,¢,,¢,,¢,, such that 


OF (6, 8, W,t,,t,,¢,) = Constant (4-10) 


4.3.2 Objective Function Jacobian 


The NAG routine EO4UCF requires the Jacobian of the objective function to be 
supplied for most efficient operation. Since the objective function has been defined as 
a constant, the Jacobian, which is the first partial derivative of the objective function 


in each independent variable is identically equal to zero in all six variables. | 


°This simple but profound idea was first suggested by Dr. Nikiforos Papadakis of the MIT 
Ocean Engineering Design Laboratory. He has conducted extensive research in optimization 
methods. 
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So the Jacobian of the objective function, J/g, may be defined as 


dOF 
0b 
OOF 
08 
OOF 
oy 
UO) oya | = 
ot, 
sor| | 
dt, 
OOF 
ot, 


(4-11) 





1S 62 SO > 


—— 








Therefore the Jacobian for the constant objective function is the zero vector. 


4.3.3 Constraint Function 


The constraint function consists of inequality constraints on the oriented 
distance function given as 

d(R,,Q;) = n, + (R;- Q) (4-12) 
If q; is the projection of a transformed measured point r; defined by (4-8), then a 

new oriented distance function after a transformation operation will be given by 
d(r,q,) = 1, + (t;—q) (4-13) 
At each iteration a new transformed measured point r;, and a new corresponding 
minimum distance projection q;, are determined for each measured point R; existing 
before the transformation. This procedure produces a new oriented distance d (r;,q;) 
for each measured point at each iteration, using the most recent transformed point as a 


Starting point for the computation of q,. 
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4.3.4 Constraint Function Jacobian 


The Jacobian for the constraint functions are the values of the first partial 
derivatives of the functions in each of the independent variables of the problem. The 
definition of the oriented distance function as the constraint function for this problem 
allows these derivatives to be calculated in a straightforward manner. 

Since the transformed oriented distance function is the constraint function for 
the localization problem during any particular iteration step, the determination of the 
Jacobian for this function consists of calculating a set of m first partial derivatives in 
the six variables ¢, 9, y, ¢,,¢,, and t,. The Jacobian J;, 1 Si Sm, for this problem may 


thus be defined as 


dd(r,.q,) on, * (r,—q,) an; + ([C]R;+t-q,) 
6 a6 a8 
dd(r,, q;) on, * (r,- q;) on, * ((CJR;+t- qi) 
ele) ele) 00 
dd(r,, q,) on, * (r,—4q,) on; * ((C]R, +t-4q,) 
ay ay oy 
Ji = Jade.a)| = |ane(n-q)} = fan-dciR+t-q)| 4! 
Ot, Ot, Or, 
dd(r;,q;) on; + (r;-4q,) dn; + ([C]R; +t-4q)) 
at, Or, ot, 
dd(r,, q;) on; * (r;—4q,) on; * ((C]R,; +t — q,) 
or, or, dt, 


which when the partial derivatives are expanded is equivalent to 
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an, _{ a[C] dt 04, 
30 (CJR, +t-—q,) + n; ic) od 
on; oC] ot 94, 
50 ({CJR,+t-—q,;) + n, ESS oy 76 au 
on, a[C] at st (4-15) 
—-°((CJR,+t-—q;) + n,° ESD Seo = ae 
ihe = oy oy oy oy 
i | on, a[C] ot 04; 
an -({CJR,+t-—q,) + n-( aE R, + an _ a 
on, o[C] ot oq; 
at, ({C JR, +t-4q,) + Nn, an, R, or, - | 
on; a[C] ot 04; 
ay ((CJR,+t-—q) + n, ac 12S 
It is clear that 
at at AIC] _ aC] _ aC} _ 
Geno Oe ee Ao) 1 et (4-16) 


By assumption, q; 1s an orthogonal projection of r; on P(u,v) and n, is the unit 
normal vector to the surface P(u,v) at the projection q;. Therefore the vector 
({C]R;+t-—q,;) is collinear with the unit normal vector n; and the following 


relationships exist: 








Oq; oq; oq; 
; 00 ; 30 = Qn, ow =a) (4-17) 
gq; Oq; Oq; 
e aa =o on. ot, = Qf, an = 0 (4-18) 
an, an, an, 
=— -({CJR;+t-q,) = =, ° (CJR; +t-q,) = =~: (CJR; +t-q,) = 0 (4-19) 
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on, on, 
(Gi SS eel ate ayo ((C]R;+t-q,) = a ((CJR,+t-—q,) = 0 (4-20) 


Now by using the expressions of (4-16) through (4-20) the Jacobian J, defined in 


(4-15) may be simplified to 














n° aR 
a[C] 
eae RR. 
n;° Ap, “ op 
o[C ] 
aC] n° R; 
n.-——R., 00 
i OW o[C ] 
J; = n; * A(t) — Pome at R; (4-21) 
or, i 
n; * o(t) ( 
dt, " 
n; ° (t) \ fe 
ot, 


where n,, , and n,' are the respective scalar components of the normal vector n,. The 
rotational elements of the Jacobian require calculating the first partial derivatives of 
the rotational transformation matrix [C] in each of the variables @, 8 and y. Those 


partial derivatives are given by the following three matrix equations. 


—cos 9 sind —cosOcoso 0 
= | cosywcosd—sinOsinysind -—cosysind—sinO@cosdsny O (4-22) 
sinycosd+sinOsindcosy -—sinysind+sinOcosdcosy 0 


a{C] 
db 


IC] —sin8coso sin @ sind cos 6 
a cos@sinycosd —cosOsindsiny  sinysin®@ (4-23) 
00 | . 

—cos@cosdcosy cosOsindcosy —sinOcosy 
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0 0 0 
= | -sinysind+sin@cosycosd —-sinycos-sin@sindcosy —cosycos@ (4-24) 
cos ysind + sin@cos¢sin y cosycosdé—sin@sindsiny  —cos@siny 


a{C] 
oy 


The rotational elements of the Jacobian are finally determined by multiplying 
each of the three matrices of (4-15), (4-16) and (4-17) by R; and substituting the 


results into (4-14). 


4.4 Problem Solution 


Using the constraint functions of (4-6) and (4-7) with the Jacobian of (4-14) as 
inputs to the optimaization routine, iterative solutions are calculated seeking a feasible set 
of rotations and translations to satisfy the constraints of the problem. The user provides 
the values for the constraints, € and 6, which define the limits for the oriented distance 
function for each measured point. The appropriate constraint is determined by the 
proximity or non-proximity of a measured point to the leading edge of the design surface 
blade. If the algorithm can find a feasible solution to the problem, that solution represents 
a prescription for the rotations and translations which should be performed on the set of 
measured points to localize them and satisfy the given constraints. If a feasible solution 
cannot be found, then the user will need to either relax the specified constraints or remove 
some measured points from the set that is analyzed. 

Experimental results using the constrained localization method are presented in 


Chapter 5. 
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Chapter 5 


APPLICATIONS OF LOCALIZATION 


5.1 Introduction 


The focus of this thesis is the development of reliable computational methods for 
the solution of the localization problem for application in the inspection of manufactured 
marine propellers. A theoretical basis for these methods was presented in Chapters 3 and 
4, but if the methods cannot be used for the intended application, then their theoretical 
development becomes only an academic exercise. The true value of the development can 
only be demonstrated if the methods can be used to solve real problems with real data. It 
is, therefore, important to show that the methods work using actual measured data and 
actual designs from manufactured propellers. 

At the beginning of the localization investigation it was decided that measured data 
from a real manufactured propeller was essential to the development and validation of the 
localization methods. The design description of a manufactured blade and a set of 
measured data points from that blade were needed for testing of the localization methods 
as they were developed. This need was identified at a periodic meeting of the 
PRAXITELES user’s group in October 1990.’ Since it was expected that the results of the 


localization investigation would produce an enhancement to PRAXITELES, the Applied 


"PRAXITELES is an interactive geometric modeling system for sculptured curves and 
surfaces. It has been developed in the Ocean Engineering Design Laboratory at MIT with 
funding from various U. S. government agencies [Hottel 91], [Tuohy 91]. 
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Research Laboratory (ARL) at Pennsylvania State University agreed to provide design 
and inspection data for a real propeller [Holter 90]. This design was a fan blade that had 
been designed, manufactured and inspected at ARL. 

At a subsequent PRAXITELES user's group meeting in January 1991 the 
preliminary results of the unconstrained localization method were presented. At this time 
it was suggested that a more complex test model for the localization method would be 
useful. While the fan blade from ARL was certainly a propeller, it did not have the 
complex sculptured geometry that would be typical of many marine propellers. For this 
reason, Philadelphia Naval Shipyard (PNSY) agreed to provide the investigator with data 
for another blade which would be more representative of marine propellers for 
localization testing [Koehler 91]. 

The experimental results and validation of the developed methods for unconstrained 
and constrained localization are presented in this chapter. Test results for the methods 
will be presented for both the ARL and PNSY blades. It 1s believed that the experimental 
results confirm the validity of the methods and provide a sound experimental basis for 


further development. 


5.2 Experimental Assumptions 


An essential element of the localization process is the determination of the 
minimum distance from a point to a surface and this distance calculation is related to the 
orthogonal projection of the point onto the surface. The orthogonal projections can be 
readily determined by exploiting some conditions that are assumed to exist in the normal 


inspection of manufactured propellers: 
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(1) The measured points on the manufactured surface can be expected to lie 
sufficiently near to the design surface so that the minimum distance projection point 
onto the design surface 1s unique. 

(2) The measurements are made away from the tip, root edge or trailing edge of the 

blade, so the minimum distance projection point is also an orthogonal projection of 

the measured point onto the surface. 

(3) The inspection device makes measurements by moving in an organized manner 

over the blade. The measurement might be made in linear or curved passes, but it is 

not random in its acquisition of data.’ 

If a cubic B-spline curve is fitted through a set of measured points using a least 
squares fitting routine with chord length parametrization, then the orthogonal projection 
of each fitted curve will provide a mapping of u,v parameters for each measured point in 
the parametric space of the design surface. This mapping approximates the orthogonal 
projection of each measured point having parameter ¢ with the orthogonal projection of 
the point on the least squares fitted curve having the same parameter ¢. Using these 
parameters as a Starting point, the orthogonal projection routine for a particular measured 
point will converge rapidly to the parameters in the design surface which correspond to 
the actual orthogonal projection. The method works because the measured points are 


near the design surface for the manufactured surface under consideration, and because 


* Further analysis of this assumption can be found in [Patrikalakis 90], [Kriezis 90], [Kriezis 
91], and [Wolter 85]. 


” This is not an essential assumption or a particularly strong one. The methods of 
measurement for automated inspection normally involve numerical programming of the 
inspection device that is analogous to the numerical programming of numerically controlled 
machinery. The assumption of organized paths allows for straightforward correlation of sets 
of measured data points and this correlation increases the efficiency of the localization 
process. 
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there are no discontinuities in the region of measurement. The method will not work at 
points of discontinuity or at points which are very near to the cut loci of the blade [Wolter 


85)."° 


5.3 Applied Research Laboratory Propeller 


The propeller design that was provided by ARL was received as a NURBS surface. 
The inspection data was received as x, y, z coordinates of measured points on the blade 
that was manufactured from the NURBS design description. These measurements were 
made at ARL using the Intelligent Robotic Inpection System (IRIS), which uses laser 
interferometry to obtain highly accurate measurements of surface coordinates. 1381 data 
points were received representing measurements on the pressure and suction sides of the 
blade. The design surface was a bicubic NURBS patch parameterized with 53 knots in 
the u direction and 20 knots in the v direction. The blade had a nominal radius of 12 
inches from root to tip. 

Some minor preprocessing was required to get the measured data into a form that 
was Suitable for localization with the developed methods. The first step in this 
preprocessing was to visually inspect the data received from ARL to establish a method 
for correlating the inspection data points. This was accomplished by using the 
visualization capabilities that exist in PRAXITELES. The blade and measured points from 


ARL had the appearance shown in Figure 1. 


" These unsatisfactory conditions typically exist at the trailing edge or tip of the blade, where 
the minimum distance projection of a point onto the surface of the blade is not uniquely 
defined. This fact is the basis for assumptions (1) and (2) above. 
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Figure 1. ARL Propeller Blade Showing Measured Points 
It was clear that the inspection measurements at ARL formed 10 correlated sets of 


points. These 10 sets of points were segregated into 10 groups of points through which a 


cubic B-spline curve was fit using a least squares approximation with chord length 


BAG! 


parametrization and zero internal knots.'' The orthogonal projection of the cubic B-spline 
curve onto the design surface was used as the source of u,v parameters for initialization of 
the orthogonal projection routines for individual measured points. With this initialization 
the orthogonal projection routines converge rapidly to the u,v parameters in the design 
surface for the actual orthogonal projection of each measured point. The u,v parameters 
of the projection are identical to the u,v parameters in the design surface at the point of 
minimum distance from the measured point. Using this information it is possible to 
establish a baseline file which contains the original measured points, the u,v parameters in 
the design surface at the points of minimum distance, and finally the minimum distances 
themselves. The distances at each point in the file form an initial reference condition to 


which the results of localization can be compared. 


5.3.1 Unconstrained Localization Results 


Further processing of the file containing the parametrizations of the minimum 
distance points in the design surface was performed to ensure that experimental 
assumption (2) was satisfied, i.e. to ensure that points in the interior of the blade were 
well away from the trailing edge, root edge or tip of the blade. Visual inspection in 
PRAXITELES revealed that chordwise measurements of the blade corresponded 
closely with the u parameter of the design surface, and predictably the span of the 
blade corresponded closely with the v parameter in the design surface. This condition 


made it possible to remove points from the total set of measured points which might 


'' This produces the simplest cubic B-spline, a Bezier curve. It is not necessary for the 
method to use any more knots or higher order B-spline than this. Using a Bezier curve 
provides for rapid fitting of the measured data points with sufficient accuracy to quickly find 
an orthogonal projection from each point to the design surface. 
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be too close to the tip or trailing edge of the blade. This was accomplished by 
choosing points whose projections in the design surface had u,v parameters which 
satisfied the condition (0.05 < u,v Ss 0.95). This resulted in an initial data file for 
localization containing 1261 points, which was 8.7 per cent smaller than the original 
set of data pcints. 


Root Mean Square Distance = (5-1) 





The results of the unconstrained localization and the computational time 
required for each set of points is presented in Table I.'* The translations and rotations 


for each set of points after unconstrained localization are presented in Table II. 


Table I 





Number RMS Distance RMS Distance 
of Before Localization After Localization Per Cent Computation 
Points (inches) (inches) Change Time (CPU) 
002371 001866 


RMS Distances and Computation Times 
for Unconstrained Localization of ARL Propeller 







The computations for the experimental results of this thesis were performed on a Silicon 
Graphics 4D25TG "Personal Iris" machine running at a nominal speed of 1.6 million floating 
point operations per second. 
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Table If 


(inches) (radians) 


[106 [00827 | -o0sor | ors | 0002s | -00047 | -0.0020 


Translations and Rotations 
for Unconstrained Localization of ARL Propeller 



















The transformation produced for the unconstrained localization of the ARL 
propeller produced an average reduction in root mean square distance from the 
sampled point set to the design surface of 19.8 per cent. The computational time 
required for the entire set of 1261 points was 1 hour and 7 minutes, producing a 
reduction in root mean square distance of 19.6 per cent.” 

To test the validity of the transformation obtained using the unconstrained 
localization method, the inverse of the transformation for the 21 point sample was 
applied to the original design surface. When the same points were localized to the 
transformed design surface the new transformation was the zero vector. 


PRAXITELES was then used to obtain a set of points on the transformed design 


'’ The computational time and reduction of root mean square distance for 1261 points is listed 
as a benchmark for the the performance of the unconstrained localization method on a very 
large data set. Although the results are consistent with those presented in Table I, they 
should be considered only as a measure of the time required for a very large data set. Other 
examples for very large data sets were not tested because of the long time required, and 
because data sets larger than about 100 data points did not seem to have a statistically 
significant effect upon the results of the localization process. 
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surface. When these points were localized relative to the original design surface, the 
Same transformation was produced as that using the 21 point sample of measured 
points. These results provide experimental confirmation of the Euclidean property of 
the localization transformation, as discussed earlier in Section 3.3 and shown in 


Appendix A. 


5.3.2 Constrained Localization Results 


The constrained localization problem differs from the unconstrained localization 
problem because measured points near the leading edge of a manufactured blade have 
greater influence on the localization than do points in other parts of the blade. The 
method uses the unconstrained localization as a starting point with the implicit 
assumption that global minimization of distances of measured points to the design 
surface is achieved before the start of the constrained localization. The constrained 
localization algorithm also uses the oriented distance function as a distance measure 
rather than the squared distance function. This distance function provides a highly 
accurate estimate of distance and is fundamental to the constrained localization 
method as developed in Chapter 4. 

The constrained localization algorithm was evaluated using the same datasets 
that were used for the unconstrained localization problem. The leading edge region 
of the blade was defined as the set of u,v parameters of the design surface where the 
condition (0.4 <u <0.6) existed. The non-leading edge region of the blade was 
defined as the set of u,v parameters for which this condition did not occur. This 
selection was based upon visual observation in PRAXITELES of the ARL blade. The 


leading edge of the blade was almost exactly coincident with the isoparameter line 
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u = 0.5 in the design surface of the blade. With this definition of the leading edge 
region, the measured points are uniquely mapped to either a leading edge region or a 
non-leading edge region based upon the value of the parameter u for the projection of 
each point onto the design surface at the beginning of the constrained localization 
process. Clearly the set of points in the leading edge region and the set of points in 
the non-leading edge region are complementary subsets of the universal set of 
measured points. 

Since the testing of the constrained localization method was intended primarily 
to demonstrate the viability of the method, the absolute magnitudes of the constraints 
were not considered as important as the relative magnitudes of the constraints for each 
region. Using this philosophy, the constraints were arbitrarily assigned with relative 
magnitudes having a 10 to 1 ratio. This means that the distance constraint in the 
non-leading edge region of the blade was assigned a magnitude of 10 times the 
magnitude of the distance constraint in the leading edge region. This constraint 
assignment strategy and the conditions necessary for a satisfactory constrained 
localization can be summarized in the following statements. 

(1) The constraints in the leading edge region of the blade have one-tenth the 

magnitude of the constraints in the non-leading edge region of the blade. 

(2) The condition for satisfactory constrained localization requires two 

necessary corollary conditions: 

a) All measured points that are assigned to the leading edge region of the 
blade must have post-localization minimum distances to the design surface 


which are less than or equal to the magnitude of the leading edge constraint. 
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b) All measured points that are assigned to the non-leading edge region of the 
blade must have post-localization minimum distances to the design surface 
which are less than or equal to the magnitude of the non-leading edge 
constraint. 

The procedure that was used to evaluate the constrained localization algorithm 
involved selecting a value for the leading edge and non-leading edge constraints 
below the threshold where a feasible localization solution could be obtained for a 
given set of measured data points. Observing that the fixed ratio between the 
magnitudes of the constraints was always maintained, the values of the constraints 
were incrementally increased until a feasible constrained localization solution could 
be found. This experimental procedure thereby determined a lower bound on the 
values of constraints which could produce feasible solutions to the constrained 
localization problem for the given set of measured points." The experimental results 
for the constrained localization of points near the leading edge and non-leading edge 
regions of the ARL blade are presented in separate tables for clarity. The root mean 
square distances, computed using equation (5-1), and the computation times for the 
two regions are presented in Tables III and IV.’ The corresponding translations, 


rotations and maximum distances are presented in Table V. 


'“In a manufacturing setting the magnitudes of the constraints would be specified by the 
blade designer. Those magnitudes might not have the same fixed ratio that was used in these 
experiments. 


'° The computation times shown in Tables III and IV are the same because the same sets of 
points were used for each table. They are duplicated in the two tables for easy reference. 
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Table III 


oe Leading RMS Distance RMS Distance 

Edge Before Localization} After Localization | Per Cent | Computational 
Points _Constraint_ (inches) (inches) Change a, 
Oo12in [001808 0.01014 


RMS Distances and Computation Times 
for Constrained Localization of ARL Propeller (Leading Edge) 










Table IV 


Number | Non-Leading RMS Distance RMS Distance 
of Edee Before Localization] After Localization | Per Cent | Computational 
Points {| Constraint (inches) _ (inches) Change Time co 


me rm ee 








a [colt a0mre [onan [sas | omnia 
er [orrin [ore | 00se77 | +1856 | Omin36 sec 
[106 [oarin [002657 | 003935 [135.0 | 1min 8see 


RMS Distances and Computation Times 
for Constrained Localization of ARL Propeller (Non-Leading Edge) 
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Table V 
(inches) (radians) Edge 
re ee Max Distance] Max Distance 
eee et ee 
er -c0018 | 0.0077 [0.1717 [00038] 0.0054 0.0008 | 0.01700in | 0.11643 in 


| 106 _| -0.1788 | -0.0629 | 0.6799 | 0.0034 | -0.0062] -0.0036 | 0.01700 in_| 0.12961 in 


Translations, Rotations and Maximum Distances 
for Constrained Localization of ARL Propeller 











It should be noted that the translation t, in Table V is typically much larger than 
the translations ¢, and f,. This result is due to the particular orientation of the blade 
relative to the axes of the coordinate system in which the measurements are made. 
The z-axis in this example is nearly parallel to the span of the blade, so that points at 
extreme positions of the z-axis must be near to the root edge and tip edge of the blade. 
These are the points which limit translational motion along the z-axis [Gunnarsson 
87a]. Since points very close to these edges are excluded from the data set at the 
beginning of the localization, there are necessarily relatively few to constrain motion 
in the z direction compared with the number of points that constrain motion in the x 
and y directions. 

The transformation produced for the constrained localization of the ARL 
propeller reduced the root mean square distance from the sampled point set to the 


design surface near the leading edge by an average of 35.7 per cent and increased the 
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root mean square distance from the sampled point set to the design surface away from 
the leading edge by an average of 90.7 per cent. The increase is a direct result of 
reducing the root mean square distance near the leading edge and transforming the 
point set away from a condition of global minimization. All points of each sampled 
Dont set had an absolute distance which was less than or equal to the listed constraint. 


The effects of the constrained localization on the global root mean square distances 


are presented in Table VI. 


Table VI 













Number RMS Distance RMS Distance 
of Before Localization After Localization Per Cent 
Points (inches) (inches) Change 


fea inl 


Global Localization Effects 
for Constrained Localization of ARL Propeller 


The global root mean square distance for the ARL propeller increased an 
average of 53.9 per cent for the five sampled point sets after the constrained 


localization process. 
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5.4 Philadelphia Naval Shipyard Propeller 


The propeller information that was received from Philadelphia Naval Shipyard was 
considerably more difficult to analyze than the data that was received from ARL because 
PNSY could not provide an analytic design descnption of the propeller blade in the form 
of a NURBS surface. Instead PNSY provided two sets of measurements of a 
manufactured blade. This created significant problems because the localization methods 
require a NURBS surface description of the design surface for their proper operation. 

The measured data from PNSY was obtained using the Automated Propeller Optical 
Measurement System (APOMS). It uses a laser interferometry technique similar to that 
of the IRIS system at ARL. The measurements were made on a submarine propeller 
blade that was manufactured many years ago by Philadelphia Naval Shipyard.'° The 
blade had an approximate radius of 68 inches form root to tip, and was therefore about 6 
times larger than the blade received from ARL. 

In order to properly test the methods of unconstrained and constrained localization, 
it was necessary to obtain a “design surface" from the set of measured points. This was 
accomplished by least squares fitting of cubic B-spline curves through the measured data 


points, followed by lofting a bicubic B-spline surface through the set of B-spline curves 


'© The measured data for this blade was certified by Commander, Philadelphia Naval 
Shipyard and Naval Sea Systems Command (Code 56X73) to be not of a classifiable nature 
in PNSY letter 9245, Code 266, Serial 9166002 of 6 February 1991. 
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so generated. Chord length parametrization was used for the curves with sixty knots.”’ 
The surface which was produced could be used as a ‘design surface" for localization 
testing with the implicit understanding that the surface is not the actual design surface as 
in the case of the ARL blade, nor is it necessarily the best possible surface representation 
of the measured data. The lofted surface is not identical to the surface from which the 
blade was originally manufactured; it is merely an approximation of that surface. For this 
reason the lofted surface will introduce some error into the localization process and will 
ultimately reduce the level of accuracy that can be achieved from it. 

The use of a lofted surface produced from measured data at PNSY as a "design 
surface’ arose out of necessity. It was not possible to obtain another blade with both a 
NURBS surface description of the underlying design and with inspection data of the 
manufactured blade as well. For this reason the PNSY design surface that was generated 
through lofting must be viewed as a simulated design surface. This example is not as 
good as the ARL example because the original design surface was not available for the 
corresponding inspection data.’* On the other hand, the example illustrates the usefulness 


of the localization methods quite adequately. 


'” A uniform spacing scheme was used for the distribution of internal knots and for ease in 
lofting. The curve is not the very best approximation of the data; significant error arises near 
the trailing edge region of the blade because of the sparcity of knots. Non-uniform knots 
with a higher concentration of knots near the trailing edge would improve the approximation 
for individual curves. However, this scheme would greatly increase the size and complexity 
of the lofted surface and it was not deemed necessary for the demonstration of these 
experiments. 


'* It should be noted that marine propeller design surfaces typically arise from lofting a 
surface through a set of curves which describe a hydrofoil section. The surface so produced 


is defined as the "design surface". Therefore the lofted surface of this example is a ‘design 


surface". The important distinction is that this design surface is not the same design surface 
as the one used for the manufacture of the propeller. 
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Despite the previously mentioned problems, the set of "design" points that was 
received from Philadelphia allowed a reasonably good blade surface to be generated. The 
“design” points were received in 17 sets of cylindrical measurements taken in radial cuts 
from the root to the tip of the blade. Each set of measurements contained 180 points 
distributed around the surface of the blade. The B-spline curves generated from these 
points and the bicubic NURBS surface that was lofted through them were all produced 
using PRAXITELES.” 

The other set of measured points that was received from Philadelphia had similar 
form to the 17 sets of points that were used to produce the design surface. There were 18 
radial cuts taken on the pressure and suction sides of the blade at radii extending from 0.3 
to 0.99. The points at radii above 0.90 are very close to the tip of the blade, so to ensure 
that all points used were relatively far from the blade tip, only points taken at radii from 
0.3 to 0.8 were used. These points and the lofted surface were visually inspected using 
PRAXITELES. The "design" blade surface and measured points for this blade had the 


appearance shown in Figure 2.”° 


'° The bicubic NURBS surface was generated using uniformly spaced internal knots with a 
total of 17 knots in the u direction and 56 knots in the v direction. The selection of direction 
for the u,v axes used for lofting was arbitrary. The resulting "design surface" parametric axes 
are orthogonal to the ARL design surface parametric axes by pure coincidence. 


*° As in the case of the unconstrained localization of the ARL blade the 12 "curves" of 
measured data which are visible in Figure 2 were fitted with cubic Bezier curves to facilitate 
finding the orthogonal projection of each measured point onto the design surface. 
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Figure 2. PNSY Propeller Blade Showing Measured Points 


5.4.1 Unconstrained Localization Results 


A total of 4077 measured inspection points from the manufactured blade were 


received from Philadelphia. To prevent problems with points near parametric 


boundaries of the design surface, these points were culled to a set of 3214 points by 
imposing the same condition that was used for the ARL blade. Points were chosen 
whose projections in the design surface had u,v parameters which satisfied the 


condition (0.05 < u,v <$ 0.95). This smaller set of points represented a 21 per cent 


reduction in the number of measured points. 
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Subsets of the global set of 3214 points were extracted for testing using the 
same procedures as those used for the ARL blade. Points were pseudorandomly 
selected over the pressure and suction sides to produce five sets of measured points 
for evaluation. The points were distributed over the entire surface of the blade to 
provide a good distribution of data. The root mean square distance, calculated using 
(5-1), was again used as a measure of the performance of the unconstrained 
localization algorithm. The results of the evaluation for five subsets of points are 


presented in Tables VII and VII. 


Table VII 











Number RMS Distance RMS Distance 
of Before Localization After Localization Per Cent Computation 
Points (inches) (inches) Change _Time ee 


a ea a NS ee ee 


ar 059 oats admin ae 
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RMS Distances and Computation Times 
for Unconstrained Localization of PNSY Propeller 


Table VIII 








Translation Rotation 
(inches) (radians) 


ae 
a1 | 0.0682 | -0.4959 | -0.0666 | -oo07e | 00028 | -00011 
Tai | o4se0 | -osiie | -oasse | -o0022 | 00035 | 0.0009 


Translations and Rotations 
for Unconstrained Localization of PNSY Propeller 












The unconstrained localization of the five sets of measured points reduced the 
root mean square distance from the points to the design surface by an average of 30.1 


DOMGCNt. 


5.4.2 Constrained Localization Results 


The constrained localization testing of the PNSY propeller used the same sets of 
data and essentially the same procedures as those that were used for the ARL blade. 
Points were assigned to a leading edge region or a non-leading edge region of the 
blade based upon the position of the projection of the point onto the design surface at 
the beginning of the constrained localization procedure. A procedural distinction 
between the PNSY blade and the ARL blade was the direction of the u,v parameters in 
the design surface. While the leading edge for the ARL blade was nearly coincident 
with the isoparameter line u = 0.5, the PNSY blade had the leading edge nearly 
coincident with the isoparameter line v = 0.5. This difference required a change in the 
definition of the leading edge and non-leading edge regions of the blade. The leading 
edge region of the blade was defined by the condition 0.4 Sv $0.6. The 
non-leading edge region of the blade was the region where this condition did not 
occur. Using these definitions a given set of measured points was mapped to specific 
regions of the design surface as was done with the ARL blade. 

The same assignment of the ratio of constraint magnitudes was used as in the 
ARL blade testing. The points assigned to a leading edge region had a constraint with 
magnitude equal to one-tenth of the magnitude of the constraint for points in the 


non-leading edge region of the blade. For testing of the process, the constraint for 


Ese 








each region was incrementally increased until a threshold was reached where a 
feasible solution to the constrained problem was obtained. In this way, a lower bound 
was determined for constraints which could provide feasible solutions to the problem 
for each set of points. The experimental results using this procedure are presented in 
Tables LX through XI.”! 

Table IX 


RMS Distance RMS Distance 
atten Before Localization| After Localization | Per Cent | Computational 
Points | Constraint (inches) (inches) Change Time AD 
bea LT SN 0109419 
010981 007713 
[ror [ oasin | oarza7 | 0.09010 | 199 | 1min22 sec 


RMS Distances and Computation Times 
for Constrained Localization of PNSY Propeller (Leading Edge) 












Table X 


Number | Non-Leading! RMS Distance RMS Distance 
of Edge Before Localization} After Localization | Per Cent | Computational 
Points Constraint pl inches) _ __(inches) — | Change Time a, 


ar [tin [ona onsers | visa | Omin Iaee 
a | 12m [ome [one | i718 | Omin 22500 
[sr] asin | 009593 | 020067 | +1028 | Omin 18 sec 


RMS Distances and Computation Times 
for Constrained Localization of PNSY Propeller (Non-Leading Edge) 











*! The computational times shown in Tables LX and X are the same because the tests were 
performed at the same time using the same set of points. The times are duplicated for easy 
reference. 
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Table XI 

of (inches) (radians) Edge 
Points = ia Se Pe ee ee Max Distance} Max Distance 
a1 |-03292| 02712 | 00487 [-0.0160]-0.0196]-0.0065| o1119in | 067443 in 
a1 | 0.1762 [04528 |-02197|-0.0013| 0.0080] 0.0021 | 0.14612in | 0457801 


| 01 | -0.1731 | 0.2832 | -0.1185 | -0.0122 | -0.0160 | -0.0032] 0.12856in_] 0.60269 in_| 


Translations, Rotations and Maximum Distances 
for Constrained Localization of PNSY Propeller 









The transformation produced for the constrained localization of the PNSY 
propeller reduced the root mean square distance from the sampled point set to the 


design surface near the leading edge by an average of 21.2 per cent. All points of 


aach A led DOInt ‘ had a absol = dista = which wa: lan Or e@auad Om 9{- 
listed constraint. The transformation produced for the constrained localization of the 


PNSY propeller increased the root mean square distance from the sampled point set to 
the design surface away from the leading edge by an average of 139.6 per cent. This 
increase was the result of reducing the root mean square distance near the leading 
edge and transforming the point set away from a condition of global minimization. 
The effects of the constrained localization on the global root mean square distances 


are presented in Table XII. 
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Table XII 


Number RMS Distance RMS Distance 
of Before Localization After Localization Per Cent 
Points (inches) (inches) Change 


[a] 0ia30 «YC 


[a _| 09760 fossa | 942 
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Global Localization Effects 
for Constrained Localization of PNSY Propeller 

The global root mean square distance for the ARL propeller increased an 
average of 101.4 per cent for the five sampled point sets after the constrained 
localization process. The entire set of 3214 points was not tested. 

The results of this test showed that the lower bound on the constraint generally 
increased with the number of points tested, but not in all cases. The constraint for 81 
points was higher than that for 101 points. An important consideration is that the 
results represent only a feasible solution to the problem for a given set of points. This 
solution is not intended to be a global minimization of the root mean square distance, 


nor is it likely to be unique. 


5.5 Unconstrained Localization of Multiple Surface Patches 


A current practice in the design of marine propellers is to break the blade surface up 
into several regions and to provide a separate design surface patch description for each 
individual region, rather than providing a single surface description of the entire blade. 


The blade may be broken up into separate patches for the leading edge, for the trailing 
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edge, for the pressure side, for the suction side, for hub fillets, etc. Each of these regions 
might in fact be broken into separate surfaces as well. It is, therefore, important that a 
localization method for marine propellers be able to accomodate the multiple surface 
patches that might exist in a real propeller design description. 

The unconstrained localization algorithm presented in Chapter 3 was generalized to 
accomodate multiple untrimmed NURBS patches. The two design surfaces and measured 
points of the ARL and PNSY propellers were used to test the algorithm. Each design 
surface was divided into three patches to represent the leading edge, pressure side and 
suction side regions of a multiple surface blade. These three patches were used to 
simulate the design representation of a multiple patch blade. It is assumed that all patches 
of such a blade are untrimmed NURBS patches of orders M in u, and N in vy, 


(O S u,v S 1) with at least tangent plane continuity. 


5.5.1 ARL Propeller Blade 


PRAXITELES was used to split the ARL blade into NURBS patches comprising 
approximately one-third of the original design surface each. The surface was split 
along the two isoparameter lines u = 0.33 and u = 0.66, and each patch was 
reparameterized such that (0 S u,v < 1). Points were then assigned to the three 
patches based upon the u,v parameters of the projection of each point onto the original 
design surface. To allow for problems that might arise for points very near to 
parametric boundaries, a "buffer" region of parametric values was assigned at the 


boundaries of the three patches. This buffer region was set equal to 0.02 in each 


aie 





parameter, so that each patch was reduced by this amount around its edges. Using 
these ideas the original design surface and the 1261 original measured points 


produced three untrimmed patches and point sets having the following characteristics. 


Pressure Side Patch (0.02 < u < 0.31 and 0.02 < v S 0.98) 402 points 
Leading Edge Patch (0.35 < u < 0.64 and 0.02 < v < 0.98) 394 points 
Suction Side Patch (0.68 < u < 0.98 and 0.02 < v < 0.98) 424 points 


The points assigned to the three regions of the blade were pseudorandomly 
selected to produce five sets of pseudorandom data as used in previous testing. (These 
point sets had comparable numbers of points, but were not the same point sets.) The 
root mean square distances from the sets of points to the surface patches were 
calculated using (5-1). The results of the unconstrained localization of these sets of 


data relative to the three surface patches are presented in Tables XII and XIV. 


Table XIII 


Number RMS Distance RMS Distance 
of Before Localization After Localization Per Cent Computation 
Points (inches) (inches) Change Time (CPU) 
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RMS Distances and Computation Times 
for Unconstrained Localization of Three Patches from ARL Propeller 
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Table XIV 


Rotation 
(radians) 


20 [02278 | 0067s | 02386 | o0o1t | 00080 | 0.0082 


[ar [02000 | -00801 | -0.1102 | 0.0031 | -00078 | -00035 
[er | osois | 00587 | ose | 0003s | -0.000s | -00022 





Translations and Rotations 
for Unconstrained Localization of Three Patches from ARL Propeller 
Using the three patches from the ARL blade the unconstrained localization 
algorithm produced a transformation which reduced the root mean square distance 
from the points to the surface patches by an average of 24.6 per cent. The 
unconstrained localization for a single patch which described the entire surface 
produced an average reduction in root mean square distance of 19.8 per cent using 


ifferent sets of me ints. 


5.5.2 PNSY Propeller Blade 


Testing of multiple patches from the Philadelphia blade followed a scheme 
virtually the same as that used for the ARL blade. The surface was split using 
PRAXITELES, but since the blade was oriented with the leading edge at v = 0.5 in the 


parametric space of the design surface, it was necessary to define the patches in terms 
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of the parameter v rather than u. The 3214 points of the PNSY design surface were 
assigned analogously to those in the ARL tests with the patches and point sets having 
the following characteristics. 


Pressure Side Patch (0.02 < u < 0.98 and 0.02 < v < 0.31) 857 points 
Leading Edge Patch (0.02 < u s 0.98 and 0.35 < v S$ 0.64) 1146 points 
Suction Side Patch (0.02 s u s 0.98 and 0.68 < v s 0.98) 939 points 


When subsets of data points were once again generated from these sets of 
points, the unconstrained localization algorithm was tested for the three patches of the 
PNSY blade using exactly the same procedure as that which was used for the ARL 
blade. The results of this testing are presented in Tables XV and XVI. 


Table XV 







Number RMS Distance RMS Distance 
of Before Localization After Localization Per Cent Computation 
Points (inches) (inches) Change Time (CPU) 


RMS Distances and Computation Times 
for Unconstrained Localization of Three Patches from PNSY Propeller 


-64- 





Table XVI 


(inches) (radians) 
03568 
ar] ao | asta | -aasv0_| eanie |"-a0oe | ooo 
Tar] 01358 | 06049 | 0.1327 | -00024 | 0.0035 | 0.001 


Translations and Rotations 
for Unconstrained Localization of Three Patches from PNSY Propeller 


















The localization algorithm produced a transformation which reduced the root 
mean square distance from the points to the surface patches by an average of 33.8 per 


cent. This compares with an average reduction of 30.1 per cent for the single patch 


description of the blade with different sets of measured points. 
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Chapter 6 


CONCLUSIONS AND RECOMMENDATIONS 


6.1 Summary of Results of Investigation 


This thesis has presented the theoretical development and numerical 
implementation in efficient computer codes of the unconstrained and constrained 
localization algorithms for application to the automated inspection of marine propellers. 
Experimental results using actual marine propeller designs with measured data have been 
provided to demonstrate the validity of the method. In each example, it has been shown 
that by using the rigid body transformation provided by the localization method a set of 
measured inspection points can be brought into closer agreement with a design surface. 
The unconstrained localization algorithm provides a method for minimizing the sum of 
Squares of minimum distances from the measured points to the surface. The constrained 
localization algorithm provides a method for satisfying localized constraints, so that 
points near particular regions of the design surface are brought closer to the surface than 
other points near other regions of the surface. The latter method provides special 
usefulness for the problem of localizing inspection points near the leading edge of a 


marine propeller blade. 


6.2 Projected Benefits of Investigation 


It is presumed that the products of this investigation, in the form of the 


unconstrained and constrained localization methods with associated computer codes, will 
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find direct application in the inspection of manufactured marine propellers. The methods 
are flexible in their requirements for input information, requiring only a NURBS 
description of a blade and corresponding spatial coordinates of measured points on a 
manufactured surface. These pieces of information can be readily obtained from existing 
systems in the industry. The unconstrained localization method is expected to be used in 
an investigation by Bird-Johnson Company of Walpole, Massachusetts and David Taylor 
Research Center in Carderock, Maryland under a U. S. Navy contract.” The method 
would be used to evaluate the inspection results of an actual marine propeller blade. 
Westinghouse Machinery Technology Division of Pittsburgh, Pennsylvania, the Applied 
Research Laboratory at Pennsylvania State University, Metal Working Technology 
Corporation of Johnstown, Pennsylvania and Martin Marietta Energy Systems at Oak 
Ridge National Laboratory have all expressed interest in the process for future marine 
propeller inspection applications. These applications are precisely ones which were 
anticipated during the development of the method. 

Beyond the direct application of the localization method which has already been 
implemented, it is expected that the development will prove useful in the area of better 
programming of work during propeller blade manufacturing. Specifically, it will be 
possible to better evaluate initial blade castings by confidently determining if an initial 
casting satisfies dimensional requirements. It is expected that fewer castings will be 
wasted because the manufacturing engineer will be able to determine the proper 
Orientation of a casting to "find the blade” in a casting which might have otherwise been 


rejected. It seems reasonable that the localization methods will aid in planning and 


* A preliminary users manual was produced directly from the work of this thesis to aid 
designers and manufacturing engineers in using the available tools [Jinkerson 91]. 
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execution of postcasting work as well. In particular, if a blade is placed in the localized 
position before it is attached to the hub, and if the attachment flange is made to conform 
to the localized blade, then an improvement in the conformance of the overall propeller to 
the original design can be expected. Finally, the localization method should be very 
valuable in the inspection of a blade before final acceptance. If the transformation 
returned from the localization process is smaller than specified tolerances for a specified 
set of measured points, then the blade could be considered to satisfy the inspection 


criteria for acceptance. 


6.3 Areas for Further Investigation 


Perhaps the one area of this work which shows the most need for further 
investigation is the existing problem of finding a suitable orthogonal projection of 
measured points at the boundaries of the design surface. In Chapter 5 the measured data 
points were carefully selected to ensure that their orthogonal projections would fall well 
within the parametric boundaries of the design surface. The localization methods, as they 
presently exist, will fail if this condition is not satisfied. Such a situation needs to be 
corrected. Extensive preprocessing of data is needed to ensure that the localization will 
work for points near parametric boundaries. Further development is needed to provide a 
value for the point projections at the boundaries. 

In Chapter 4 an assumption in the development of the constrained localization 
problem was that the projection of a measured point on the design surface would not 
change during the localization process. Although this assumption is almost certainly 
valid for very small rotations and translations, it constrains the flexibility and generality 


of the process. Rather than fixing the mapping of measured points to assigned regions 
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from beginning to end of the process, it would be desirable for the assignment to be 
changeable as the projection of a point may change during the localization. This problem 
is not expected to be trivial. It is expected that discontinuities will result as constraints 
change dynamically, and these may be difficult to handle with current optimization codes. 
More theoretical development will probably be necessary to address this problem. 
Alternatively, a smooth bivariate function might be employed to represent the constraints 
throughout the patch, rather than by using the piecewise constant constraints that were 
used in this work. 

The constrained localization problem in its present formulation produces a feasible 
but not necessarily unique solution. It starts from a position of presumed global 
minimization, which is hopefully achieved during the unconstrained localization process, 
and seeks a condition which will satisfy local constraints. A much more difficult problem 
is one which would gptimize the solution to the constrained localization problem. It is 
clear from the premise that the constrained problem starts from a global minimum, that 
any subsequent solution cannot be an unconstrained global minimum. Such a solution 
will however satisfy the optimality condition of constrained global minimization. A 
related problem would be to show that the unconstrained localization does in fact produce 
a global minimum. 

Another area of investigation that would be worthwhile involves the application of 
the localization methods, constrained and unconstrained to the problem of trimmed 
NURBS patches. Untrimmed patches have been addressed in this work, a trimmed 
patches will be more difficult. As a minimum, a satisfactory solution to the problem of 
finding the orthogonal projection of a measured point at the boundaries of a design 


surface must be found. Solution of the trimmed patch problem will, however, greatly 
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increase the generality of the localization methods. 

The accuracy and efficiency of the constrained localization algorithm presented in 
Chapter 4 was improved by providing symbolic Jacobians for the constraint function and 
objective function. It is likely that the unconstrained localization algorithm could be 
improved by supplying the Jacobian or higher order derivatives to the minimization 
routine as well. The Jacobian for the squared distance function is derived in Appendix B. 
An obvious next step in the future development of the localization process is to 
implement it in the unconstrained localization algortihm. 

A final area that is worthy of investigation involves the application of statistical 
theory to the selection of measured points for evaluation and to the results produced from 
the methods. While the simple root mean square distance measure is an appropriate one 
for macroscopic evaluation of the results of the method, it would be very interesting to 
apply some statistical measures to the results in order to evaluate the specific effects of 
points in various regions of the blade. Furthermore, the selection of numbers of measured 
points to be used in the evaluations of the examples of Chapter 5 was consistent but 
arbitrary. Statistical experiments should be performed to determine the appropriate size 


of point sets for a given set of measured data. 
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Appendix A 


DEMONSTRATION OF THE ORTHOGONALITY 
OF THE ROTATIONAL TRANSFORMATION MATRIX 


A 3 by 3 matrix [A] is defined to be orthogonal if it can be shown that 


[A][A)’ = [A] [A] = 4, (A-1) 


where /; is the identity matrix of dimension 3. 

It was asserted in Chapter 3 that the rotational transformation matrix presented there 
was orthogonal. It will be shown through a direct, albeit somewhat tedious, application of 
matrix algebra that equation (A-1) 1s satisfied and that the rotational transformation matrix is, 
in fact, orthogonal. This fact has significance in the context of this work because it shows 
that an inverse transformation of the rotational transformation matrix, identically equal in this 
case to the transpose of the matrix, would return any point in space that was operated upon 
by the matrix to its original position. The matrix then can represent a valid geometrical 
transformation. 


The rotational transformation matrix [C] 1s defined in (3-10) as 


cos 8cos —cos 8sind sin 8 
[C] = | cosysind+sinOsinycosh cosycosd—sinOsindsiny —sinycos® (A-2) 
sinysind—sinO@cosdcosy sinycosd+sinOsindcosy cos8cosy 
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[C}’ is then the transpose of [C], defined as 


cosO@cosd cosysind+sinOsinycosd sinywsind—sinOcosdcos y 
[C]’ = | -cos@sind cosycosd—sin@sindsiny sin ycosd+sin@sindcos y (2) 
sin§9 —sinycos 9 cos 8cos y 


Let the following assignments be made: 


a=cosy c =cos8 e=cosd 
b=sny d=sin9 f=sing (A-4) 


Using these assignments, the two matrices [C] and [C]’ become 


ce —cf d 
[C] = |af+dbe ae-dfb —-bc (A-5) 
bf-—dea be+dfa ca 
ce af+dbe bf—dea 
[(C! = | -cf ae-—dfb be+dfa (A-6) 
d —bec ca 
Now let [A] = [C][C]’ and [B] = [C}’[C] 
Proceeding with the matrix multiplications, 
A,, = (c’e’) +(c’f’) +(a’) 
A, = (acef + bcde’) + (-acef + bcdf’) + (—bcd) 
A,, = (beef — acde’) + (—bcef — acdf’) + (acd) 
A,, = (acef + bcde*) + (—acef + bcdf” )+ (—bcd) (A-7) 


A, = (a’f’ + 2abdef + b’d*e’) + (a’e’ — 2abdef + b*d’f’) + (b’c’) 

A,, = (abf —a°def + b*def — abd’e’) + (abe* + adef — b’def — abd’ f’) + (-abc’) 
A,, = (bcef - acde’) + (—bcef — acdf” )+ (acd) 

A, = (abf® + b’def — a°def — abd?e”) + (abe? — b°def +.a°def — aba?f*) + (-abc?) 
A,, = (bf — 2abdef +.a°d’e’) +(b’e* + 2abdef +a°d’f’) +(a°c’) 


Jon 





B,, =(c’e’) +(a°f’ + 2abdef + b’d’e*) + (b’f’ — 2abdef +a°d’e’) 

B,, = (-c’ef) + (a’ef —abdf + abde* — b*d’ef) + (bef +abdf —abde’ —a’d’ ef) 
B,, = (cde) + (-abcf — b’cde) + (abcf —a°cde) 

B,, = (—c’ef) + (a’ef + abde* — abdf’ — d*b*ef) + (b’ef — abde* + abdf —a’d’ef) 
B,,=(c'f)+(a°e’ — 2abdef + b*d’f’) +(b’e* + 2abdef +a°d’f’) 

B,, = (—cdf) + (-abce + b’cdf) + (abce +a*cdf) 

B,, = (cde) + (-abcf - b*cde) + (abcf —a*cde) 

B,, = (—cdf) + (-abce + b*cdf) + (abce +a°cdf) 

B,, =(d°)+(b'c’)+(a’c’) 


(A-8) 


Using the trigonometric identity cos’ a+ sin’ a = 1, then 
a+b =c’+d’ =e’ +f =) (A-9) 


The equations in (A-7) and (A-8) can be simplified using these identities. 


A, =Cetcped@acl(e’+fPyac'+d’=1 
Aj, = acef + bcde’ — acef + bcdf’ — bcd = bcde’ + bcdf’ — bcd = bcd(e’ +f’ -1) =0 
A,, = (bcef — acde’) + (-bcef — acdf’) + (acd) = -acde’ — acdf’ + acd 
= acd(-e’-f'+1)=0 
A,, = (acef +bcde’) + (-acef + bcdf’) + (-bcd) = bcde’ + bcdf’ - bcd 
= bcd(e’+f’-1)=0 
A,, = (a°f" + 2abdef + b’d’e*) + (a’e’ — abdef + bd’ f’) + (b’c’) 
sa ftb'de tae’ +b dftb'c’ 
= arle’+f)+b'd'(e*+f’) +b°c=a*t+b'd'+b'c’= a’+b*(c’+d’) =a’+b’=1 
A,, = (abf’ — a’def +b’def - abd’e’) + (abe’ +. a’def — b’def - abd’ f’) + (-abc’) 
=ab(f-d’e*+e?-d’f)—c’ 
=ab{1-d°(e’+f’)—c’] = ab[1 -(c’ +d’)] =0 
A,,= (bcef - acde’) + (-bcef — acdf’) + (acd) = —acde’ - acdf’ + acd 
= acd{1-(e*+f")] =0 
A, = (abf’ + b’def — a’def — abd’e’*) + (abe’ — b’def + a’def — abd’ f’) + (-abc’) 
= abf’ - abd’e’ + abe’ — abd’ f’ — abc’ 
=ab(f'-d’e’+e°-d’f'-c’) =ab{l —d'(e* + f")-c’] =ab[1 -(d’+c’)] =0 
A,, = (b’f? - 2abdef + ad’e’) + (b’e* + Zabdef +.a’d’f’) + (a’c’) 
=b’fta'd'e’+be’+a'd'f tVcrader+f)ta'd'(e?+f)t+a'c’ 
=b+a'(c’+d’)=a’ +b’ =1 


(A-10) 


ie 








B,, =(c’e’)+(a’f’ + Zabdef + b’d’e’) + (b’f' — Zabdef +.a’d’e’) 
=Ce ta feb eth f+a'd'e’ 
=e that Dt death yaclet+fedeae(ctdytfeetfeal 

B,, = (-c’of) + (a’ef — abdf’ + abde* — b’d’ ef) + (bef + abdf’ — abde’ — ad’ ef) 
=-c’ef +a°ef —b'd’ef + bef —a°d’ef =-c’ ef + ef (a° +b’) -d’ef(a* +b’) 
= ef{1—(c’ +d’)] =0 

B,, = (cde) + (-abcf — b’cde) + (abcf —a°cde) = cde —b’cde —a’cde 
= cde{1—-(a’?+b*]}=0 

B,, = (-c’ef) + (a’ef + abde’ — abdf’ — d*b’ ef) + (b’ef + abde’ + abdf? — ad’) 
=-c ef +a’e —d’b’e +b’ eo —a'd'e =ef(a’+b’)—-d'e(a’+b’)-c’ef (A-11) 
= ef[1-(c’+d’)}=0 

B,, = (c’f°) + (a’e’ — 2abdef + b’d’f) + (b’e* + Yabdef +a’d’f’) 
=CPt@etbh da ftberadfeacfr+e(a’ +b +a f(a’ +b’) 
setf(ctd =e +f =l 

B,, = (-cdf) + (-abce + b’cdf) + (abce +a°cdf) = —cdf + b’cdf +a’cdf 
= -—cdf + b’cdf +a’cdf = cdf(a’ +b’) -cdf =0 

B,,= (cde) + (-abcf — b*cde) + (abcf —a’cde) = cde — b’cde —a’cde 
= cde —b’cde —a’cde =cde[1—(a’ +b’)] =0 

B,, = (—cdf) + (-abce + b’ cdf) + (abce + a’cdf) = —cdf + b’cdf +a’ cdf 
= cdf[(a’+b’)- 1] =0 

B,=a(@)+(Po)t(Vo ad’ +b'c +a’ ad +c'(a’ +b) =c? +b =1 


From (A-10) and (A-11), 


A,,=8,,=1 A,.=8,,=0 A,,=8,,= 
A,,=B,,=0 A,.=B,= A,, =B,, = 


(A-12) 
A,,=B,,=0 A;=B,,=0 A,,=B,;= 
It is clear then that 
SIP Sk See SG eh ae (A-13) 


The matrix [C | is therefore orthogonal. 


QED 
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An alternate proof of the orthogonality of the rotational transformation matrix, i.e. that 


equation (A-1) is satisfied can be based upon the fact that 


[C] = [EJ [FIG] 


(A-14) 


where [E], [F] and [G] represent the three matrices on the right hand side of (3-10). It is easy 


to show that these matrices are orthogonal. 


1 0 0 cos8 QO sin@\fcos@ -sinod 0 
(C] = |0 cosy -siny 0 | 0 sind cosd QO 
QO sinyw cosy /\-sin@ O- cos@ 0 0 ] 

l 0 0 1 0 0 


[EJ(EY =|0 cosy -siny|/0 cosy siny|=4J, 
O siny cosy /\O -siny cosy 


cos® Q sin8@\{cos8 O -—sin@ 


[FI[FY =| 0 1 + O Oo 1 0 |=4, 
—-sin8 QO cos@/\sin8 QO cos8 
cosd -—sind O)\{ cosd- sind QO 

[G][Gy' =| sind coso O|]| -sind coso 0] =1, 


0 0 ] 0 0 Nt 


Then using matrix properties from linear algebra: 


(EJIFIIG)’ (EMFIUG) = (GY (FY (EY (E] (FG) 
(GY (FY (EY (E](FJIG) = (GY (FY (EY (E) (FIG) = (GY (FY [FIG] 
(GY (FY (FJ(G] = (GY (FY (F)(G] = (GY [G] = J, 


The commuted expression is shown similarly. 


QED 


-75- 


(A-15) 


(A-16) 


(A-17) 


(A-18) 


(A-19) 


(A-20) 


(A-21) 





Although this proof is straightforward, and vastly simpler than the proof outlined in 
(A-2) through (A-13), the former proof is useful because it provides an independent 


verification of the symbolic correctness of (A-2). 
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Appendix B 


DERIVATION OF JACOBIAN 
FOR SQUARED DISTANCE FUNCTION 


The Jacobian of the squared distance function which is used as the objective function in 
the unconstrained localization problem of Chapter 3 is estimated in the NAG routine E04JAF 
by difference quotients. Because it 1s expected that another minimization routine could make 
productive use of a symbolic Jacobian, the Jacobian for the squared distance function is 
presented here. Such a routine could be the routine EO4UCF used in Chapter 4, with 
application to the unconstrained localization problem for enhanced precision and efficiency. 
The derivation is analogous to that of the oriented distance function which was used as the 
constraint function in the constrained localization problem of Chapter 4. 

Recalling the notation of Chapters 3 and 4, assume again a parametric design surface, 
P(u,v), a set of m measured points R,, and the set of nearest points or projections, Q,. The 


minimum distance from a measured point R; to the design surface P(u,v) is defined as 
d(R;,Q;) = | R;-Q;| = min,,,d [R;, Pu, v)] (B-1) 


A new set of transformed points r,, can be obtained from the original set of measured 
points R;, by the following operation: 

r. = (CJR, +t (B-2) 

If r; is a transformed point given by equation (B-2), and q; is the projection of the 

transformed point onto the design surface, then the transformed minimum squared 


distance from r; to q,, {d(r;,q,;)]’, may be defined as 


oe a 





(d(r,,q,)) = (r,- 4) ° (r;-@) (B-3) 


The Jacobian for the squared distance function is the set of first partial derivatives of 
the squared distance function in the six independent variables of the problem. Since the 
transformed minimum squared distance function is the objective function for the localization 
problem during any particular iteration step, the determination of the Jacobian for this 
function consists of calculating a set of m first partial derivatives in the six variables 9, 9, y, 


t,, t,, and t,, The Jacobian J;, 1 Si Sm, for this problem may thus be defined as 


ld(r:,a)T A(r; — 4) + (r;- q;) 
do ob 
dld(r,,q,)F o(r;— q;) * (t;-q;) 
00 00 
o[d(r,,q,) o(r; — q;) * (Fr; - q;) 
y = oy = oy (B-4) 
! o[d(r;, q,)]’ O(r; — q;) * (t;— gq) 
Or, Ot, 
d{d(r;,q;) o(r; — q;) * (tr; - q;) 
dt, Oty 
ld (r;.4)F Ari =.) * (i= 4) 
or, dt, 


which when the partial derivatives are expanded is equivalent to 
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Cle , 24) 
Z [Ae + dd as | (r; q;) 
Cl, , AM) 
2 (ae 00 i * 90 3) pee 
(Ay més a), 4) 
ay ay 
Jy = (B-5) 
2| a =R, + ed 
| Ot, Ot, 74 
0 ot oq, 
2 ( ~ a 4) tw 
F) ot oq; 
2 ean eae 4) 
It is clear that 
ot _ ot _ ot AC] _ AC] _ AC) _ 
ob 60608) OW o7 ane or, Oty ot, 0] (B-6) 


As developed in Chapter 4, by assumption, q; is an orthogonal projection of r; on 
P(u,v) and n, is the unit normal vector to the surface P(u,v) at the projection q;. Therefore 


the vector (r; —q,) 1s collinear with the unit normal vector n, and the following relationships 


exist: 
oq; oq; oq; 
00 (r;-q,) = 30 ad) oy (r;-q) = 0 (B-7) 
0q; - oq; = = 0 
ag (6 = ae x, (r; - q)) leed:) = (B-8) 


By using the expressions of (B-6) through (B-8) the Jacobian J; of (B-5) may be 


simplified to 


-79- 








BY us (r, — q;) 
Cle sg _ 
2 ae R; (r; q;) 
o[C] 
2 ——R.°(r.-q. 
= ae (r;—q,) 
; = ot t Zz (B-9) 
Or, i q; 
ot 
2 or, : (r;-q;) 
2 Eva (r;—4q,;) 
and this may be finally simplified to 
2 2 R; * (r;- q;) 2 0 R; * (r; - q,) 
2 76 R; + (F; — q;) 2 By R; * (r; — q,) 
J; = 2 ay Bw) a ay Rn (B-10) 
Ze (Gah) 20 =.) 
2e\ + (r,;-q;) 2(r, - 4) 


Lee diag.” JE e200 —4.) Oey 
In this expression e,’, e, and e,' are unit vectors along the axes of the coordinate system which 
are parallel to the components of the translation vector t. The three differences r,'-g,', r,'-q, 
and r,'-q,' are the scalar components of the difference vector r,-q;. The rotational elements of 
the Jacobian require calculating the first partial derivatives of the rotational transformation 
matrix [C] in each of the variables @, 8 and yw. Those partial derivatives are given by the 


following three matrix equations. 
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IC} —cos 8singd — cos 8cos ° (B-11) 
oy) = | cosycosod—sinOsinysind -—cosysind—sinOcosdsiny 0 
sinycosd+sinOsindcosy —sinysind+sinOcosdcosy 0 
IC] —sin@coso sin 8 sind cos 8 (B-12) 
ier cosOsinycosd —cos@sindsiny — sinysin®8 
—cos8cosdcosy  cosO9sindcosy -—sinOcosy 
IC] 0 0 0 (B-13) 
ak = —sinysind+sinO@cosycosd -sinwcosd—sn@sindcosy —cosycos8 
¥ cos Wsind +sin®@ cos $sin y cosywcosd—sinO@sindsiny —cosOsny 


The rotational elements of the Jacobian are then determined by multiplying each of the 
three matrices of (B-10), (B-11) and (B-12) by R; and substituting the results into (B-10). 

The Jacobian of (B-7) was not used in this investigation because the routine that was 
implemented for the unconstrained localization problem could not use gradients supplied by 
the user. It is expected that the use of this Jacobian would greatly improve the efficiency and 
accuracy of the optimization process, and it is clear from this development that the Hessian 


and other higher order partial derivatives could be readily derived. 
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Appendix C 


DETERMINATION OF PARAMETERS IN DESIGN SURFACE 
FOR HIGH ACCURACY 
MINIMUM DISTANCE CALCULATION 


The very accurate determination of minimum distance from a point to a parametric 
Surface that is used in the oriented distance function of Chapter 4 uses the u and v parameters 
in the design surface of the orthogonal projection of the given point onto the design surface. 
Therefore the problem is to find as accurately as possible the pair of u, v parameters which 
correspond to the projection of a given measured point onto the design surface. 

In Chapter 3 a modified Newton algorithm (NAG routine EO4KCF) was used to 
determine the minimum distance from a measured point to the design surface for use in the 
squared distance function. This method provides a fast and reasonably accurate 
determination of the minimum distance by finding a minimum of the squared distance of the 
measured point from an arbitrary point P(u.v) with parameters u,v in the vicinity of a starting 
point approximation (uy, ¥,). Given the inherent inaccuracies of any minimization routine in 
floating point arithmetic, the solution will be only approximate.” 

An improvement in the accuracy of the calculation of minimum distance from a 
measured point to the design surface can be achieved by using the orthogonality of the 
projection to determine the values of the parameters uw and v in the design surface. If the 


values of u and v from the minimization routine are used as a starting point, then the 


= Using the NAG routine EO4KCF in 16 digit floating point arithmetic, only 7 digits of 
precision in the calculation of the paramters u and v for the minimum distance calculation 
can typically be obtained. 
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orthogonality condition, which by an assumption in section 3.2 must exist at the projection of 
the measured point onto the design surface, provides a simple method for obtaining a very 
accurate value for the minimum distance. The development of the method is similar to that 


presented in Appendix C of [Kriezis 90]. 


Using the notation of Chapters 3 and 4, consider again the minimum distance d(R;, Q,), 


from a measured point R,, to the design surface P(u,v). Since Q, is defined as the projection 


of R; onto the design surface P(u,v), orthogonality necessarily requires that 


(R;-Q,) > 0,P(u,v) = 0 (C-1) 
and 
(R;-Q;)* 0,P(u,v) = 0 (C-2) 
where 
Q; = P(u,v) (€2) 


Given the necessary and sufficient conditions for orthogonality of (C-1) and (C-2), the 
determination of the u and v parameters in the design surface P(u,v) which correspond to the 
projection point Q; can be accomplished by finding the zeros in u and v which satisfy the 


conditions of (B-1) and (B-2). So consider two functions F ,(u,v) and F (u,v), defined as 


F, (u,v) = (R; — Q;) . 0, P(u,v) = (R; — P(u, v)) i 0, P(u, v) (C-4) 


and 


F,(u,v) = (R;-Q;) + 0,P(u,v) = (R;- Plu,v)) + 0,P(u, v) (C-5) 


ce 





The zeros in u and vy for functions F (u,v) and F,(u,v) are found using the NAG routine 
COSPBF which utilizes another modified Newton method to find the zeros of multivariable 
functions with the gradients in each variable supplied by the user. 


The gradients for F,(u,v) and F (u,v) can be expressed as 


0,F, = -d,P(u,v)*0,P(u,v) + (R;-P(u,v)) + 0,,P(u,v) (C-6) 
0,F, = -0,P(u,v) + 0,P(u,v) + (R,-P(u,v)) + 0,,P(u,v) (C-7) 
0,F, = -0,P(u,v) + 0,P(u,v) + (R;- P(u,v)) * 0,,P(u,v) (C-8) 
0,F, = -0,P(u,v) + 0,P(u,v) + (R;-—P(u,v)) + 0,,P(u,v) (C-9) 


There is obviously a time penalty associated with the use of this method rather than the 
simple minimization routine. The improvement in accuracy may justify the use of this 
method when accuracy is more important than computational speed. In particular, this 
method was employed in the oriented distance function computation used for the constrained 


localization problem of Chapter 4. 
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